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ABSTRACT 



The contribution of the Source Cosmic Rays (SCRs), confined in Supernova 
Remnants, to the diffuse high energy 7-ray emission above 1 GeV from the 
Galactic disk is studied. 7-rays produced by the SCRs have a much harder 
spectrum compared with those generated by the Galactic Cosmic Rays which 
occupy a much larger residence volume uniformly. SCRs contribute less than 
10% at GeV energies and become dominant at 7-ray energies above 100 GeV. 
The contributions from 7r°-decay and Inverse Compton 7-rays have comparable 
magnitude and spectral shape, whereas the Bremsstrahlung component is 
negligible. At TeV energies the contribution from SCRs increases the expected 
diffuse 7-ray flux almost by an order of magnitude. It is shown that for the inner 
Galaxy the discrepancy between the observed diffuse intensity and previous 
model predictions at energies above a few GeV can be attributed to the SCR 
contribution. 



Subject headings: gamma-rays - background radiation - cosmic rays - supernova 
remnants 
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1. Introduction 

Observations of the diffuse Galactic 7-ray emission give information about the 
Galactic Cosmic Rays (GCRs), the interstellar gas and diffuse photon fields, and about the 
interactions between them. The observational results obtained with the Energetic Gamma 
Ray Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory can be 
described fairly well by a suitable model for the diffuse interstellar gas, Cosmic Ray (CR) 
and photon distributions (e.g. Hunter et al. f997a; Hunter et al. 1997b). 

However, above 1 GeV the observed average diffuse 7-ray intensity in the inner Galaxy, 
300° < I < 60°, \b\ < 10°, exceeds the model prediction significantly. There are at least two 
possible explanations for this discrepancy (e.g. Hunter et al. 1997b; Weekes et al. 1997). 
The high-energy 7-ray excess may indicate that the GCR spectrum observed in the local 
neighborhood is not representative of the diffuse CR population in the Galactic disk; a 
harder average diffuse proton spectrum is required to explain the 7-ray excess if it is due to 
7r°-decay. An unresolved distribution of CR sources is the other possibility. 

The physical picture which we consider in this paper corresponds to the second 
possibility. The idea that CRs, after leaving their sources, could in principle produce 
7-rays in ambient dense clouds with a harder spectrum than those produced by the average 
GCRs, was proposed by Aharonian and Atoyan (1996). Another proposed possibility, 
also invoking transport effects, is the local hardening of the CR energy spectrum in the 
direction perpendicular to the Galactic disk above strong CR sources, especially above the 
inner region of the Galaxy, due to faster CR convection in a faster Galactic Wind (Volk 
1999). This contributes a principally observable harder than average 7-ray component in 
such regions. In contrast, we consider here the accelerating particles inside their sources, 
where they are much more strongly scattered than in the ISM, neglecting the contributions 
invoked by Aharonian and Atoyan, and by Volk (see also below). 
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We assume that SNRs are the dominant sources of the Galactic CRs. On this premise 
we find that CRs, accelerated and confined in SNRs, give an important contribution to 
the high-energy 7-ray emission from the Galactic disk. Since the CR energy spectrum 
inside SNRs is much harder than on average in the Galaxy — the average spectrum 
being softened by rigidity-dependent escape from the Galaxy in the diffusion region above 
the disk — the relative SNR contribution increases with energy and becomes in fact 
dominant at 7-ray energies e 7 ^100 GeV. It may substantially increase the diffuse TeV 
7-ray emission from the Galactic disk so as to constitute a significant and spatially variable 
observational background which must also be taken into account in the search for spatially 
extended Galactic CR sources in this energy region. A physically analogous problem is the 
contribution of CR electrons in SNRs to the radio synchrotron spectrum of normal galaxies, 
without an Active Galactic Nucleus. It has been recently discussed by Lisenfeld & Volk 
(1999). 

In this paper we shall investigate pion-decay 7-ray emission from CR nuclei as well 
as Inverse Compton (IC) radiation and Bremsstrahlung due to CR electrons. It will be 
shown that the average IC 7-ray background from SNRs is comparable in magnitude and 
spectral form to the pion-decay background at high energies, whereas the corresponding 
Bremsstrahlung component is negligible. 



2. Gamma-ray luminosity of old SNRs 

The majority of the GCRs, at least up to kinetic energies e ~ 10 14 eV, is presumably 
accelerated in SNRs. According to modern theory a significant part of the hydrodynamic 
Supernova (SN) explosion energy E S n ~ 10 51 erg is converted into CRs already in the 
early Sedov phase of the evolution, due to diffusive shock acceleration (e.g. Berezhko et al. 
1996; Berezhko & Volk 1997). Later on, the CR energy content and the high-energy 7-ray 
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production slowly decrease with time. This is at least true as long as the progenitor star 
is not so massive as to have a strong wind which significantly modifies the circumstellar 
medium (Berezhko & Volk, 2000). The total number of SNRs N$n = ^snTsn is an 
increasing function of their assumed life time T$n, i-e. the time until which they can confine 
the accelerated particles; here v S n is the Galactic SN rate. Therefore we conclude that 
the population of the oldest SNRs dominates the total 7-ray luminosity of the ensemble of 
Galactic SNRs. Thus we consider only old SNRs which nevertheless are still strong enough 
to confine most of the CRs produced during the prior evolutionary stages. 

We then have the situation that the CRs in the Galaxy are represented by two basically 
different populations. The first one consists of the ordinary GCRs and presumably occupies 
a large Galactic residence volume quasi- uniformly. This residence volume exceeds that 
occupied by the CR sources by far (e.g. Ptuskin et al. 1997; for an earlier review, see 
Berezinsky et al. 1990). The second CR population, which we call Source Cosmic Rays 
(SCRs), is represented by shock accelerated CRs that are still confined in the localized 
SNRs. During the initial, active period of SNR evolution of about 10 5 yr when the SN 
shock is relatively strong, the volume occupied by the accelerated CRs practically coincides 
with the shock volume. In later stages the shock becomes weak and CRs begin to leave the 
SNR acceleration region. After some period of time T SN the escaping SCRs become very 
well mixed with the "sea" of GCRs. We shall assume that the transitory period during 
which SCRs are transformed into GCRs is much shorter than the preceding confinement 
time Tsn- The most important factor for CR confinement is the shock strength. Even in 
the phase where radiative cooling would formally become important, this remains true since 
CRs prevent cooling compression due to their pressure. Thus the assumed mean life time 
Tsn~ 10 5 yr is determined by the shock dynamics more than by anything else. 

Since the 7-ray production due to GCRs is quite well studied (e.g. Hunter et al. 
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1997b; Mori 1997), it is primarily important to find the relative contribution of the SCR 
population. 



2.1. Gamma rays from 7r°-decay 

The production rate of 7r°-decay 7-rays from inelastic CR - gas collisions, primarily 
p-p collisions, may be written in the form (Drury et al. 1994) 

Q 7 (e) = Z y a pp cN g n(e), (1) 

where N g is the local gas number density, a pp is the inelastic p-p cross-section, Z 7 is the 
so-called spectrum-weighted moment of the inelastic cross-section, n(e)de is the CR spatial 
number density of CRs in the kinetic energy interval de, and c is the speed of light. Thus 
we have to primarily calculate n(e) for the two CR populations. 

The quasi-uniform GCR population in the gas disk is assumed to have roughly a power 
law spectrum in the relativistic range 

mc z \mc z / 

The total number n^ CR of relativistic GCRs with e > mc 2 , per unit volume, can be 
expressed in terms of the CR energy density eacR' 

OCR _ (jGCR ~ 2)e G CR , v 

\1GCR - 1)^C 2 

where m is the proton mass. For simplicity we restrict our consideration here to the proton 
component which is energetically dominant in both the GCR and the SCR populations. 

In contrast to the GCR population, the SCRs are confined inside a discrete number 
Nsn of SNRs. These are assumed to be predominantly located in the Galactic gas disk, 
of volume V g . Spatially averaged over the Galactic disk volume their 7- ray production 
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rate is determined by an expression analogous to eq. (1), where instead of n(e) one should 
substitute the SCR distribution 

n S cn{e) = N S c R (e)N SN /V g , (4) 

with NscR(c)de being the overall (i.e. integrated over the SNR volume) SCR number in 
the energy interval de, and should use an appropriate mean local interstellar medium (ISM) 
number density Ng CR into which the SNe explode. 

Since the CRs produced inside SNRs have also a power-law spectrum N S cr oc e~ lsCR 
in the relativistic range, n S cR^) can be expressed in the same forms (2) and(3), with 

e S CR = N SN SEsN/Vg, (5) 

and putting ^scr > 2. Here 5 is the fraction of the SN explosion energy E$n converted into 
SCRs. 

However, according to the prediction from nonlinear kinetic theory (Berezhko et al. 
1996), diffusive shock acceleration produces an extremely hard spectrum of SCRs at the 
early Sedov phase which is characterized by a power law index ^scr = 2. In this case we 
have 

SCR _ e scR (6) 

mc z ln(e max /mc 2 ) 

where e max is the maximum SCR energy. 

For the ratio R = Q^ CR /Q^ CR of the 7-ray production rates due to SCRs and GCRs, 
we have 

Z* cr N sn 5Esn 

{tl) ZGCR( lGCR - 2 ) He max /mc*)V g e GCR 

x c(A0 , (7) 

\mc z / 

where ( is the ratio Ng CR /N g 3CR , and Ng 3011 denotes the average gas density in the disk. 
In fact we assume that the gas and the CRs are distributed uniformly inside each SNR 
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which is approximately true for the old SNRs which we consider here. The parameter ( 
describes a possible spatial correlation between SN occurrence and local ISM density. If on 
average Supernovae explode in a denser than average medium in the Galactic disk, then 
C > 1, whereas ( < 1 in the opposite case. 

The main mass of the ISM M g = 4 x 10 9 M Q is contained in a Galactic disk region 
of a thickness of about 240 pc, corresponding to the thickness of the HI gas (Dickey and 
Lockman 1990) which has the volume V 9 = 2.5 x 10 66 cm 3 referred to before. Here we take 
a disk radius of about 10 kpc which implies an average gas density N^ CR = 2 cm~ 3 . Taking 
the relativistic part of the GCR spectrum to be characterized by tccR — 10~ 12 erg/cm 3 , 
and 'Jgcr = 2.75 which results in Z^ CR /Z^ CR = 10 (Drury et al. 1994), we obtain for the 
standard set of SN parameters E S n = 10 51 erg, v S n = 1/30 yr _1 : 



in addition using the rather moderate parameter values 5 = 0.1 and e max = 10 5 mc 2 to 
characterize CR acceleration inside SNRs (e.g. Berezhko et al. 1996). One can see from 
this expression that for T$n ~ 10 5 yr the 7 ray production due to SCRs becomes dominant 
already at energies e 7 >10 GeV. 

Note that the quantity 5E SN N SN j (V g eGCR) represents the ratio of currently existing 
total SCR energy and GCR energy inside V g . For the above set of parameters it is about 
0.1. Despite the fact that the SCRs represent only a relatively small fraction of the total 
CR energy content even in the disk, they may dominate the 7-ray production at sufficiently 
high energies due to their much harder spectrum. 

It is clear that the quantity -R(e 7 ) determines the average ratio 
(dN^ CR /de 1 ) j \dN? f CR /<ie 7 ) of 7-ray spectra produced in any region of the disk, by SCRs 
and GCRs, respectively. Therefore the total 7-ray spectrum measured from an arbitrary 
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Galactic disk volume is expected to be 

dN dN GCR 

ff = ^ 4 + ^)], (9) 

where the additional factor 0.4 is introduced to approximately take into account the 
contribution of GCR electron component to the diffuse 7-ray emission at GeV energies 
(e.g. Hunter et al. 1997b). In Fig.l we present the expected differential flux of 7-rays from 
the inner Galaxy, calculated for £ = 1 and T SN = 10 5 yr, with the spectrum dN GCR /de 7 
taken from the paper by Hunter et al. (1997b), and extended into the region e 7 > 30 GeV 
according to the law e~ 2 - 75 . One can see that after inclusion of the 7-rays produced by 
SCRs, the calculated flux even exceeds the EGRET flux for e 7 >20 GeV. This suggests that 
expression (7) overestimates the 7-ray production by SCRs. 

It is possible that the overall source spectral index Jscr is somewhat larger than 2 due 
to very late accumulation of only low-energy particles. It can also not be excluded that 
the confinement time T S n of the SCRs depends on energy. Due to their high mobility, the 
highest energy particles may leave the vicinity of parent SNR earlier and also more rapidly. 
This process of SCR escape into the ISM starts for the most energetic particles already 
at the early Sedov phase of SNR evolution (e.g. Berezhko et al. 1996; Berezhko & Volk 
1997). Therefore at time T SN , when the main part of SCRs are released from the SNR, 
their spectrum may be somewhat steeper than a Jscr = 2 spectrum. 

Due to the importance of the problem we derive the relation between risen and ugcr 
in a different form. It leads to the same results if SNRs are the GCR source. We start from 
the usual leaky box balance equation 

^GCi?(e) N SC r , ir .x 
= T/ v SN , (10) 

T c Vc 

where V c (e) is the energy- dependent residence volume occupied by GCRs that reach the gas 
disk during their mean residence time r c in V c . In the case of an extended Galactic halo 
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due to a Galactic wind driven by the GCRs themselves, and for energies much larger than a 
few GeV, V c can be much greater than V g . In fact V c (e) oc e a55 in such a selfconsistent halo 
model (Ptuskin at al. 1997). 

Note that the leaky box model deals with a CR distribution ncci?( e ) averaged over 
the residence volume V c . Therefore it can only be applied to the GCRs which can be 
assumed to be almost uniformly distributed in the residence volume. It is not valid for 
the SCRs, because their behavior is determined not only by large-scale transport but also 
by other physical factors which, for example, lead to their acceleration. Technically the 
volume Vscr, occupied by the SCRs, should be excluded from the residence volume V c and 
the SCRs appear in the balance equation (10) for the GCRs only in the form of a source 
term N S crVsn/V c as a CR population released from the source region into the ISM after 
some unspecified evolutionary period T S n- Therefore eq.(10) does not depend upon T S n- 
However, the effects produced by the SCRs confined inside the ensemble of simultaneously 
existing SNRs, for example the additional 7-ray production, essentially depends on Tsn, 
since it directly determines the total number N$n of simultaneously existing SNRs. 

Using eq. (4) we can write 

nscR _ VcTsn _ Tsn 
n G CR V g r c Tg 

The GCR residence time in the disk volume, r g = r c V g /V c , can be derived from the 
measured grammage x, which is the mean mass of Interstellar matter traversed by GCRs of 
speed v in the course of their random walk in the Galaxy: 

r 9 = (12) 
The measured grammage at high energies e > eo = 4.4 GeV is 

x — 14 (v/c) (e/4.4 GeV)~ a60 g/cm 2 ; for e < e , x = lAv/c g/cm 2 (Engelman et al. 1990). 
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Therefore the residence time in the gas disk can be written in the form 

r g = T (e/6 )-», (13) 

where r = 4.6 x 10 6 yr, e = 4.4 GeV, \x = 0.6, and e > e . This experimentally inferred 
value for /i closely agrees with the theoretical result of Ptuskin et al. (1997). 

We note that at relativistic energies e > mc 2 , according to the initial balance eq.(10), 
the GCR spectrum uqcr oc e~ lGCR and the overall SCR spectrum N$cr oc e~~ isCR should 
be connected by the relation 

1SCR = 1GCR ~ H- (14) 

For \i = 0.60 the source should produce a SCR spectrum with Jscr = 2.15, while for the 
case t scr = 2 one would need /i — 0.75. 

At the same time the SCR distribution 
^5Ci?(e) oc e~ 7 scfl ; averaged over the gas disk, can have a different shape compared to 
Nscr(^) in the case of an energy dependent SNR confinement time T S n(^), according to 
eq. (11). It can be steeper, r y' SCR = 7gcr — f 1 + P, if T$n oc e _/3 is a decreasing function of 
energy. It would mean that the highest energy SCRs leave the parent SNR faster than the 
lower energy SCRs. 

Eq.(ll) leads to a simple expression for the 7-ray production ratio 

*(*r) = ^ip^' ( 15 ) 
independently of u$n- Substituting the residence time in the form (13), and taking 
7scr = 2.15, which leads to Z^ CR /Z° CR = 7.5 (Drury et al. 1994) in eq.(15), we obtain for 
e 7 > 4.4 GeV: 

We shall only consider 7-ray energies e 7 > 4.4 GeV. One can see here again that the 
SCR contribution is determined by the value of the confinement time inside SNRs, T S n- 
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Unfortunately there is no detailed description of when and how CRs, accelerated in SNR, 
are released into the ISM. Nevertheless one can give some constraints on this process. 

According to the standard theory, the expanding SNR shock produces a power law CR 
spectrum up to the maximum energy (Berezhko 1996; Berezhko et al. 1996; Berezhko & 
V61k 1997) 

e m oc R S V S , (17) 

which is determined by the radius R s and speed V s of the shock. The CRs with the highest 
energy e max are produced at the very beginning of the Sedov phase t ~ t when the product 
R S V S has its maximum RqV , where 



are the sweep-up time, sweep-up radius and initial mean ejecta speed respectively; M e j 
denotes the ejecta mass, and p g = N g m the ISM density. Subsequently, the product R S V S 
decreases with time as t -1 / 5 and the SNR shock produces CRs with progressively lower 
cutoff energy e m (t) < e max = £m(to)- During that phase those CRs that were previously 
produced with energies e m < e < e max now propagate outward diffusively without significant 
influence of the SNR shock. If their expansion is still governed by the Bohm diffusion 
coefficient as during their acceleration, the expansion rate is only slightly higher than the 
SNR expansion rate and these particles should be considered as confined inside the source 
(i.e. the SNR). 

The opposite extreme case corresponds to the assumption that particles with energies 
e > e m do no more produce a high level of turbulence. Let us consider this pessimistic 
scenario in terms of confinement here. In this case the propagation of these very high 
energy particles is governed by the mean Galactic diffusion coefficient which is very much 
larger than the Bohm diffusion coefficient. In this situation particles with energy e should 
be considered as released from the source at the moment t > to when e m (t) drops below e. 
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Since the particles with maximum energy are produced at t ~ t , one can write 

T SN (e) = t (— ) 5 . (19) 

Due to this strong dependence it is clear that T SN (e) will still deviate from the overall gas 
dynamic life time T s ^ only for large energies e near e max . All particles with 

will remain confined until T s °^. 

Since the majority of Galactic SNe are core collapse SNe from stars with masses 
exceeding about 8M , we shall use M e j = 10 M & for purposes of estimate. Except 
for progenitor masses exceeding 15 M , the progenitors have only a weak stellar wind. 
Therefore the assumption of a uniform circumstellar medium remains reasonable for the 
average properties of the SNR population. 

For the main fraction of CRs the confinement in SNRs terminates when the SNR 
shock becomes weak and produces CRs with a very steep spectrum that cannot anymore 
excite a high level of turbulence near the shock front. If we take M = 4 as a critical 
Mach number, the corresponding SNR age will be t — (M /M) 5/3 = 1.5 x 10 3 t , where 
M = V /cso is the initial shock Mach number and c S o is the ISM sound speed. For an ISM 
with number density N g = 1 cm -3 , c,so — 4 km/sec, and then to ~ 1.5 x 10 3 yr which gives 
T SN ~ 2 x 10 6 yr for M ej = 10 M and E SN = 10 51 erg. This estimate shows that the SNR 
shock remains rather strong during a very long period of time. 

Another physical factor which can restrict the SCR confinement is the radiative cooling 
of the postshock gas. Approximately it becomes important when the postshock temperature 
drops below ~ 10 6 K, or when the postshock sound speed drops below C52 ~ 100 km/s. 
In the case of a strong shock, with Mach number M ^> 1, the postshock sound speed is 
determined by the shock speed Cs2 ~ V^V s /3. During the Sedov phase the shock speed 
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decreases with time according to the law V s = 0.4 x Vq (£/£q) 3 ^ 5 - Therefore the shock 



SN and ISM parameters gives t c ~ 6 x 10 4 yr. Since gas clumping as a result of cooling may 
lead to effective SCR leakage from the SNR, a value of the confinement time T$n = 10 5 yr 
is reasonable for N g = 1 cm -3 . (Fig. 1). 



In Fig.l we present a calculated 7-ray spectrum based on the above expression (16) 

with 



tmax = 10 5 GeV, and e 7 = O.le, which is roughly valid for the hadronic considered 7-ray 
production process considered. At GeV energies in this case SCRs contribute about 10% of 
the total 7-ray flux. Due to their hard spectrum this contribution progressively increases 
with energy and becomes dominant at e 7 >100 GeV. It increases the expected TeV 7-ray 
flux by about a factor of ten. Note that the actual SCR contribution from the inner part of 
the Galaxy is somewhat higher than the above estimate due to the larger SNR concentration 
in this region. 

As one can see from Fig. 1, the SCR contribution in the case Jscr — 2 , using eqs. (8) 
and (9), is for all energies larger than that for Jscr — 2.15, using eq. (16). This difference is 
related to the different SCR acceleration efficiencies. In the first case it is characterized by 
the parameter 5 = 0.1 which, according to eq. (5), directly determines the 7-ray production 
rate. In the second case the SCR acceleration efficiency 5 = [J^S X eiVsc#(e) de]/ ' Esn is not 
contained in the final expression (15), but one can derive it easily from the balance equation 
(10) which gives: 



speed drops to the value V c = 100 km/s at the age t c = t (0.4 V /V c ) 5 ^ 3 . The above set of 



EDITOR: PLACE FIGURE [I] HERE. 



T SN = mm{10 5 , 10 3 (e/e ma ,r 5 } yr, 



(20) 




(21) 
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To determine 5 we use the (demodulated) GCR distribution 

n GCR {e) = 8.1 x 1(T 10 

X (l^V + llv) -- 3 ^V)-' (22) 
(Ryan et al. 1972; Perko 1987). 

Substituting the values of T g (e), V g , Esn, and i>sn, we obtain 5 ~ 0.05 for the 
relativistic part of the GCR spectrum. The required acceleration efficiency is two times 
lower compared with the case ^scr — 2 due to the essentially steeper spectrum. Note, that 
in both cases the SCR spectra with ^scr — 2, 8 — 0.1 and ^scr = 2.15, 5 = 0.05 contain 
about the same number of relativistic CRs. In the first case the required GCR residence 
time is r g oc e _a75 , whereas the second case with r g oc e~ 06 is close to the experiment 
(Engelman et al. 1990). Therefore we believe that the dashed line in Fig.l represents the 
most reliable estimate for the expected diffuse 7r°-decay 7-ray emission, especially at high 
energies e 7 > 100 GeV. 

Note also that the acceleration efficiency required by eq. (21), using an assumed 
Galactic SN rate v SN = 1/30 yr _1 and a mean SN explosion energy E SN = 10 51 erg, 
is considerably lower than that predicted by shock acceleration theory, which gives 
5 = 0.2 0.5 (Berezhko et al. 1996). Yet, in contrast to the acceleration models which 
determine 5 from the injection rate and the nonlinear acceleration theory selfconsistently, eq. 
(21) determines only the product 5E SN u SN from observed quantities. The observationally 
inferred SN explosion energies E S n can be at least by a factor 2 smaller than 10 51 
erg, and from comparisons with galaxies similar to our own v$n could vary between 
1/30-r- 1/100 yr _1 . Therefore the empirical value of 5 from eq. (21) can vary between 0.05 
and 1/3. Nevertheless, the theoretically determined efficiencies appear systematically too 
high. 

As a possible solution for this discrepancy one might assume that, just before being 
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released, the SCRs lose an important part of their energy by adiabatic expansion so that 
the SCRs' energy content inside a SNR is higher than the energy contained in the released 
spectrum Nscr(c)- However, there is little dynamical basis for such an assumption. Much 
more likely is that the very efficient CR acceleration inside SNRs predicted by the nonlinear 
kinetic theory, assuming spherical symmetry, takes place in reality only at some fraction 
of the SN shock surface, because suprathermal positive ion injection into the acceleration 
process on the highly oblique part of the shock can be significantly suppressed (Bennet & 
Ellison 1995; Malkov & Volk 1995). In this case the acceleration efficiency, calculated for a 
spherical SNR shock, should be reduced by a factor of a few. 

The actual SNR distribution can in fact be rather nonuniform within the disk volume, 
contrary to what we have assumed implicitly up to now. In this case the estimated value 
of i2(e 7 ), which describes the relative SCR contribution to the diffuse 7-ray emission, 
should be corrected by a factor Ng N /N S N, where Ng N is the expected number of SNRs in 
the observed region and N$n represents this number in the case of uniformly distributed 
SNRs. It is clear that the expected value of i?(e 7 ) is almost independent of the actual SNR 
distribution if the observed region is an essential part of the whole disk volume V g . 

2.2. Inverse Compton and Bremsstrahlung gamma-rays from SCR electrons 

Electrons, once being injected into the diffusive shock acceleration process, will be 
as efficiently accelerated in SNRs as are the protons. Even though there exist theoretical 
concepts (e.g. Levinson 1994; Galeev et al. 1995; McClements et al. 1997; Bykov & 
Uvarov 1999), electron injection is not completely understood. However, there is no 
doubt that electrons undergo continuous acceleration during SNR evolution. The spectral 
shape of accelerated electrons Ng CR (e) inside SNRs deduced from radio-observations on 
average agrees with what is expected from shock acceleration. Since relativistic electrons 
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with energies e > 1 GeV are dynamically indistinguishable from protons, their source 
spectrum Ng CR (e) = K ep N SC R(e) can differ from that of the protons N S cr(^) only by some 
energy independent factor K ep that is determined by the injection process. High energy 
electrons produce 7-ray emission due to IC scattering, especially on the Cosmic Microwave 
Background (CMB) and by Bremsstrahlung on the interstellar gas. We shall first consider 
the IC contribution here. 



2.2.1. Inverse Compton contribution 

In an approximate form, valid if the generating electron energy distribution is smoothly 
varying, like in the case of a power law considered here, the IC 7-ray emissivity Q 7 c (e 7 ) 
can be written as (e.g. Longair 1981, Berezinsky et al. 1990) 

Qf (e 7 ) = a T cN ph n% CR {e e )^ (23) 



where 

e e = m e c 2 Y / 3e 7 /(4e pA ) (24) 

is the energy of electrons which produce an IC photon with mean energy e 7 , 
or = 6.65 x 1CT 25 cm 2 denotes the Thomson cross section, e p h = 6.7 x 1CT 4 eV 
and N p h = 400 cm -3 are the mean energy and number density of the CMB photons, 
respectively. Finally n e SCR = Ng CR Ng N /V g denotes the average spatial electron SCR 
number density. 

Thus the ratio of the IC to the 7r°-decay 7-ray production rate reads as 
Q? _mo B , Acm-W e, ^ N* SN 



QSCR 



where we have used Jscr = 2 and a pp = 4 x 10~ 26 cm -2 ; Ng N denotes the number of SNRs 
which contribute to the IC emission at energy e 7 from CR electron sources. 
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From this expression it appears as if the IC 7-ray contribution would be dominant at 
TeV-energies if K ep is as large as 1CT 2 and if the mean gas number density inside SNRs 
is about N^ CR = 1 cm 3 . However, the radiative cooling time r e (e e ) of electrons, which 
produce 7-rays with energy e 7 = (e e /17.1 TeV) 2 TeV, 

Te = 7.3 x 10' {-£-) [—) yr, (26) 

reaches the above assumed overall SNR confinement time T S n = 10 5 yr for 7-ray energies 
e 7 < e* = 5.4(5/10 /zG)~ 4 GeV. Here B is the magnetic field strength inside the source, 
whose typical value inside SNRs is about 10 /iG. Therefore, taking into account the obvious 
relation Ng N /N SN = r e /T SN , on average the relative IC contribution of the electron 
component of SCRs in TeV 7-ray can be written as 

It is independent of the 7-ray energy for e 7 > e*, and only given by eq. (25) with 
Ng N /N S N = 1 for e 7 < e*. This consideration shows that for the parameters assumed, and 
for K ep of the order of 10~ 2 , we have an IC contribution to the average 7-ray background 
which is comparable to the hadronic background for all energies above a few GeV. 



2.2.2. Bremsstrahlung contribution 

At high energies, e e ,e 7 ^> m e c 2 , we have for the Bremsstrahlung 7-ray emissivity Q 7 r 

Qf (e 7 ) = 2 / de e -^cN s g CR n% CR {e e ), (28) 

where e m i n is the minimum electron energy necessary to produce a Bremsstrahlung 7-ray 
of energy e 7 , the factor 2 takes into account the contributions of electron-electron and 
electron-proton collisions, and where the differential electron-proton Bremsstrahlung 
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cross-section is given by (e.g. Berezinsky et al. 1990) 
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Here a 1/137 and r = 2.818 x 10~ 13 cm denote the fine structure constant and the 
classical electron radius, respectively. 

For our chosen value ^scr = 2, the integral for Q 7 r (e 7 ) can be calculated in closed 
form. In the limit ln(e 7 /m e c 2 ) 1, of interest here, it reduces to the asymptotic form 
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Thus, finally, we obtain 
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This small ratio implies that Bremsstrahlung 7-rays play no role for the average 7-ray 
background above GeV energies, if K ep <0.1, taking into account, that N§ n /Nsn is always 
smaller than 1. 



3. Discussion 

We note that of order ten SNRs of age younger than 10 5 yr can on average lie within 
a 1 degree field of view of a detector directed towards the Galactic Center. Therefore 
a moderate fluctuation of the measured 7-ray intensity is expected due to variations of 
the actual number of SNRs within such a detector's field of view. At the same time, for 
directions perpendicular to the Galactic plane, the chance to observe the contribution of 
SCR 7-ray emission is quite negligible. A question is then how one might best study the 
nonuniformities of this background experimentally. Clearly this is an investigation of its 
own. Therefore we would like to restrict ourselves to a few comments here. Due to the 
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spectral form of the background its graininess is most pronounced at high 7-ray energies. 
This is even more true due to the fact that for individual sources with a very low magnetic 
field the IC emission could be much stronger than the 7-ray emission due to hadronic 
interactions; in addition the IC emission has a harder spectrum. This suggests the use 
of imaging atmospheric Cherenkov telescopes. Their resolution in angle and energy is as 
good or better than that of other ground-based detectors. However, the study of extended 
sources is not an easy task with imaging telescopes which have a very limited field of view, 
even employing the stereoscopic method. For low brightness extended sources a satellite 
instrument like the future GLAST detector is well suited since it does not have to deal with 
the charged CR background due to its use of an anticoincidence shield. On the other hand, 
the statistics achievable with a small area space detector gets very low above a few tens of 
GeV. Thus one should probably attempt such a study with both types of instruments due 
to their complementary properties. 

A different question concerns the limitations of our approach due to its concentration 
on SNRs as the sources of the GCRs. In fact the considerations in this paper can be applied 
to any alternative class of dominant GCR sources. The most important aspect, which leads 
to the dominance of SCRs in high-energy 7-ray production, is that the GCR sources should 
generate SCRs with a spectrum that is significantly harder than the GCR spectrum. Eq. 
(15) is valid for an arbitrary class of CR sources if we substitute some other value of the 
SCR confinement time T5 instead of Tsn, since the grammage x is an experimentally fixed 
quantity. Let us then assume that some class of compact CR sources produces an energy 
Eq in the form of CRs with spectrum Nscr e~ lsCR with average frequency ug, and 
let us further assume that this spectrum remains unchanged inside the source regions for 
some period of time T s after which these CRs are released into the ISM as the GCRs. It 
is obvious that due to the general energy requirement the production rate v s E c should be 
about the same as v SN bE SN . The SCR energy E c , deposited in some initial volume V, will 
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produce a dynamically significant disturbance in the background ISM if we assume that 
the initial SCR energy density E c /V is much greater than the thermal ISM energy density 
which in turn is of order eocR- This will inevitably lead to the confinement of these SCRs 
inside an expanding, disturbed volume for some period of time Ts before the SCRs will be 
released. It is difficult to give a general relation between E c and T s and there may exist 
only lower bounds on Ts, given Ec- Therefore, we cannot exclude speculative source classes 
with many weak but long-lived sources. The opposite case of many weak and short-lived 
sources is excluded to the extent that the present explanation of the hard 7-ray spectrum 
by the contribution of the SCRs is unique. 



4. Summary 

Our considerations suggest that the SCRs can provide an essential contribution to the 
high-energy Galactic 7-ray flux. According to our estimates, depending on the parameters, 
the SCR contribution is less than 10% of the GCR contribution at GeV energies and it 
dominates at energies greater than 100 GeV due to its essentially harder spectrum. This 
conclusion is confirmed by calculations performed for the case when SNRs are the main 
source of GCRs. At TeV energies the SCRs increase the expected 7-ray flux from the 
Galactic disk by almost an order of magnitude. 

The single physical parameter which determines the SCR contribution due to hadronic 
interactions is the SCR confinement time Tsn- As far as the 7-ray emission due to 7To-decay 
is concerned, the above conclusions are valid for Tsn ~ 10 5 yr. Since this SCR contribution 
is proportional to T S n, it would be negligible at TeV energies if T SN < 10 4 yr. A SNR age 
of 10 4 yr typically corresponds to the intermediate Sedov phase, when the SNR shock is 
still quite strong. Therefore it seems to be quite improbable that the GCRs are replenished 
from SNRs at such an early phase. For the IC contribution even a ten times shorter source 
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life time would be sufficient at TeV energies. In fact, for the TeV IC emission the relevant 
time scale is the life time r e (e e ) of parent SCR electrons due to their synchrotron losses, 
which is indeed about 10 4 yr. For decreasing 7-ray energies r e (e e ) increases beyond 10 4 yrs, 
and therefore a source life time of this magnitude would become a limiting factor. 

Our conclusions remain valid for alternative classes of possible GCR sources with 
comparable overall energy release and comparable individual confinement times. We note 
that this contribution of the dominant GCR sources necessarily exists. As we argue, it may 
be sufficient by itself to explain the observed 7-ray excess, at least in the inner Galaxy 
where it is well documented, without a need to invoke additional particle populations (e.g. 
Pohl & Schlickeiser 1991). 
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Fig. I. — The average diffuse 7-ray spectrum of inner Galaxy (300° < / < 60°, \b\ < 10°). 
The full (dashed) line represents the calculation with SCR spectral index ^scr — 2 
{iscR — 2.15), and the dash-dotted line corresponds to the n° decay 7-ray spectrum 
produced by GCRs (Hunter et al. 1997b). EGRET data are also taken from the review 
paper by Hunter et al. (1997b). 



